Complex Molecules in the LI 157 Molecular Outflow 



Hector G. Arce 

Department of Astronomy, Yale University, P.O. Box 208101, New Haven, CT 06520-8101 

hector . arceOyale . edu 

o : 

CN ' Joaqum Santiago- Garcia 

^ \ Observatorio Astronomico Nacional, Alfonso XII 3, E-28014 Madrid, Spain 

\^ . j .santiagoOoan.es 

^ Jes K. J0rgensen 

^ : 

Oh' Argelander-Institut fiir Astronomie, University of Bonn, Auf dem Hiigel 71, 53121 Bonn, 

O ' Germany 

to ! 

, jestSastro.uni-bonn.de 

y—i ■ Mario Tafalla 

> : 

^ ■ Observatorio Astronomico Nacional, Alfonso XII 3, E-28014 Madrid, Spain 

vn ■ 

CN I m.tafalla@oan.es 

in ; 

^ . and 



OO 
O 



Rafael Bachiller 



^ ' Observatorio Astronomico Nacional, Alfonso XII 3, E-28014 Madrid, Spain 



r . bachillerOoan . es 



ABSTRACT 



We report the detection of complex organic molecules in the young proto- 
stellar outflow L1157. We identify lines from HCOOCH3, CH3CN, HCOOH and 
C2H5OH at the position of the Bl shock in the blueshifted lobe, making it the 
first time that complex species have been detected towards a molecular outflow 
powered by a young low-mass protostar. The time scales associated with the 
warm outflow gas (< 2 x 10^ yr) are too short for the complex molecules to have 
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formed in the gas phase after the shock-induced sputtering of the grain man- 
tles. It is more hkely that the complex species formed in the surface of grains 
and were then ejected from the grain mantles by the shock. The formation of 
complex molecules in the grains of low-mass star forming regions must be rela- 
tively efficient, and our results show the importance of considering the impact 
of outflows when studying complex molecules around protostars. The relative 
abundance with respect to methanol of most of the detected complex molecules 
is similar to that of hot cores and molecular clouds in the galactic center region, 
which suggests that the mantle composition of the dust in the LI 157 dark cloud 
is similar to dust in those regions. 

Subject headings: ISM: molecules — stars: formation — ISM: jets and outflows 
— ISM: individual (L1157) — astrochemistry 

1. Introduction 

Shocks from protostellar outflows heat and compress the surrounding medium thereby 
triggering different processes (e.g., grain disruption, ice gain mantle sublimation, etc.) that 
can release molecules trapped in the grains into the gas phase. In this way, outflows con- 
tribute to the chemical enrichment of the gaseous environment surrounding young stars. This 
is clearly demonstrated by multi-molecular line observations reflecting enhanced molecular 
abundances either in shocks associated with the outflows or in the outflow cavity walls in the 
protostellar envelopes (e.g., Garay et al. 1998; Bachiller et al. 2001; Arce & Sargent 2006; 
J0rgensen et al. 2007). In most cases these molecules have similar velocity and structure as 
the CO lines that trace the outflows, and it is clear that the outflow is responsible for their 
overabundance. 

The limited outflow chemistry studies show shock-triggered overabundance of parent 
(or "flrst generation") species — molecules that are released directly into the gas phase from 
the icy dust mantles, like H2CO and CH3OH — as well as other simple molecules presumably 
formed in the warm gas (e.g., SO, HCO"*"). Chemical models indicate that larger organic 
molecules (made of 7 atoms or more) can form on grain surfaces (e.g., Garrod & Hcrbst 2006) 
and should be present in the gas phase along with the other simpler molecules observed near 
outflow shocks — as suggested by Chandler et al. (2005) and Remijan & HoUis (2006). Yet, up 
to date, conclusive evidence for shock-triggered overabundances in molecular outflows from 
low-mass stars (M < 2M0) exists only for species with six atoms or less. Here we report 
the flrst detection of HGOOH, GH3GN, HGOOGH3, G2H5OH in the the L1157 molecular 
outflow. 
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The LI 157 dark cloud (at a distance of about 440 pc) harbors an embedded low- 
mass stellar object, L1157-mm (IRAS 20386+6751), that powers a young molecular outflow 
(Umemoto et al. 1992). The powerful L1157 outflow is comprised of several shocks seen as 
bright knots in near and mid infrared images (Davis & Eisloeffel 1995; Cabrit et al. 1998; 
Looney et al. 2007). Millimeter observations of many molecular emission lines show that 
L1157 is a chemically rich outflow exhibiting an overabundance of SiO, CH3OH, H2CO, CN, 
HCN, SO, and a number of other species by a factor of a few tens to more than 10^, de- 
pending on the molecule (Bachiller & Perez Guitierrez 1997, hereafter BP97). In the blue 
(southern) lobe, the emission from these species is brightest towards two to three clumps 
(Bachiller et al. 2001), and high- resolution interferometer observations reveal the clumps are 
associated with shell-like structures created by two jet-driven bow shocks (Zhang et al. 1995; 
Gueth et al. 1996). These results indicate that the observed chemical enhancements are 
induced by the outflow shocks. 



2. Observations and Results 



Observations of the hnes shown in Table 1 were obtained using the IRAM 30 m telescope 
in Pico Veleta, Spain, in June 2007. Our observations concentrated in the brightest CO clump 
in the blue lobe of the L1157 outflow (i.e., the Bl clump), located at an offset of Aa = 
22.5", A6 = -64.5" from the L1157-mm protostar at a(J2000) = 20''39"06".2, 5(J2000) 
= 68°02'16" (Bachiller et al. 2001). The Bl clump is far enough from the source and the 
telescope beam is sufficiently small (see below) that any molecular line emission associated 
with the gas immediately surrounding L1157-mm (i.e., within 10^ AU) does not contaminate 
our observations. Two different spectral configurations were used in order to include as 
many fine transitions as possible of HCOOH, C2H5CN, CH3CN, and HCOOCH3 (A and 
E) — complex molecules that have been detected in the immediate surroundings of deeply 
embedded low-mass protostars (Cazaux et al. 2003; Bottinelli et al. 2007). Both spectral 
configurations included four receivers, each simultaneously connected to a 256 channel filter 
bank with a spectral resolution of 1 MHz and a unit of an autocorrelator, that provided 
simultaneous observations at four different frequencies ranging from 90.1 to 257.6 GHz. The 
autocorrelator was set to provide spectral resolutions of 40 and 320 kHz and bandwidths 
between 40 and 320 MHz. All observations were obtained in wobbler switching mode with a 
90" throw — large enough to avoid any possible "contaminating" emission from the rest of the 
blueshifted lobe or the immediate surroundings of LI 157- mm. Verification of the telescope 
pointing approximately every two hours, indicated that the pointing accuracy was between 
2 and 4". A total integration time of four to ten hours, and average system temperatures 
between 120 K and 830 K, depending on the frequency, resulted in spectra with rms of 0.002 
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to 0.01 K. The telescope beamwidth (FWHM) at 90.2, 98.6, 135.7, 226.7, and 257.5 GHz is 
27", 25", 18", 11", and 10", respectively. Spectral data presented here are in units of main 
beam brightness temperature, (see Rohlfs & Wilson 2000). 

We detect, at a 3 cr level or higher, five HCOOCH3-A transitions, four HCOOCH3-E 
transitions, three CH3CN transitions, two HCOOH transitions and one C2H5OH transition 
(see Table 1). These were identified using the CDMS and JPL catalogs (Miiller et al. 2001; 
Pickett et al. 1998). We are confident that the emission arises from the LI 157 outfiow's blue 
lobe as all lines are blueshifted by 1.5 to 4.0 km s~^ with respect to the LI 157 cloud central 
velocity {Vlsr = 2.7 km s~^) and have a velocity width between 2 and 7 km as expected 
from the velocity distribution of other species associated with the LI 157 outflow, observed 
towards Bl by BP97. Figure 1 shows an example of our detect lines. We used the rotational 
diagram method to estimate the rotational temperature (T^ot) ^-nd column densities of the 
species for which we detect two or more transitions (see, e.g., Bisschop et al. 2007). We only 
detect one transition of C2H5OH, so we estimate its column density assuming Trot ~ 25 K 
(similar to that of HCOOCH3) and follow Requena- Torres et al. (2006). For all lines we 
assume the emission is optically thin and the gas is in local thermodynamic equilibrium. 
In Table 2 we show our results with a correction for beam dilution, assuming a source size 
of 10.5". We base our assumption of the source size on the high-resolution interferometer 
CH3OH observations of the L1157 blueshfited lobe by Benedettini et al. (2007) which show 
that Bl is mostly made of two subclumps. Taken together, the emission from these two 
subclumps extends about 10.5" (estimated from the geometric mean of the minor and major 
axis of the emission, see Figure 2 of Benedettini et al. 2007). Whether the complex molecules 
are formed as a by product of methanol in the warm gas or are ejected from the grain mantle 
with the methanol (see below), we expect the emission from the observed complex molecules 
to be coincident with that of CH3OH (see also, Liu et al. 2001). The CH3CN emission is not 
corrected for source size as the beam size for all three transitions of this molecule is smaller 
than 10.5". 

Our estimates indicate that most of the detected molecules have Trot between 10 and 
30 K. This is lower than the estimated gas temperature of Bl (T ~ 80 — 300 K, see below), 
as well as significantly lower than the temperature needed to evaporate off most of the 
species present in icy grain mantles (e.g., van Dishoeck & Blake 1998). Even so, similar low 
temperatures have also been derived for the same complex molecules in the gas immediately 
surrounding low- mass protostars (Cazaux et al. 2003; Bottinelli et al. 2004; 2007). It is 
possible that the low estimates of Trot are due to our poorly constrained rotational diagrams, 
as we only detect two to four transitions (mostly at low energies) for each molecule. In 
addition, non-LTE effects may infiuence the derived temperature from the rotational diagram 
(Goldsmith & Langer 1999). Shock-heated molecular gas takes only a few hundred years to 
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cool from 100 K to 20 K (Bergin et al. 1998), and it is plausible that some of the detected 
molecules come from gas that has cooled down to the estimated low temperatures. Further 
observations are needed to distinguish with certainty which of the scenarios suggested above 
affect our estimates. 

We calculate the abundance of each molecule {X — N/NH2) using our estimate of the 
molecule's column density (A^) and the estimate of H2 column density {Nh^) towards Bl 
obtained by BP97. The latter was obtained from single-dish observations of high-velocity 
CO outflow emission and assuming a standard CO/H2 ratio of 10^^. Uncertainties on the 
size of Bl, the CO/H2 ratio and the opacity of the CO line introduce uncertainties on 
our estimate of the absolute value of the abundance. However, these uncertainties are less 
important when comparing the relative abundance between the different molecules. 

3. Discussion 

Our results clearly indicate that complex organic molecules can be found in outflows from 
low-mass stars. The interesting question that follows is, how are these complex molecules 
formed? The two main competing scenarios propose that complex organic molecules mainly 
form either in the warm gas phase or on grain surfaces. The gas phase models were originally 
conceived to explain the chemical richness of hot cores. In these objects the radiation 
from a massive protostar is believed to heat the dense {n ~ 10^ — 10^ cm~^) inner core to 
temperatures higher than 100 K where the icy grain mantles evaporate injecting molecules 
hke H2O, H2CO, CH3OH into the gas-phase (e.g., Kurtz et al. 2000). The gas phase models 
indicate that subsequent chemical ion-molecule reactions in the warm gas can result in the 
formation of other, more complex, molecules, such as HCOOCH3 and C2H5OH (the so 
called daughter or "second generation" species) (Charnley et al. 1992; Caselli et al. 1993). 
Alternatively, in the grain surface model, the observed complex organic molecules are formed 
on the grain surface and are later released into the gas phase (with H2CO and CH3OH) when 
the grain mantles evaporate (e.g., Hasegawa et al. 1992; Hasegawa & Herbst 1993). 

In both models, a heat source is necessary to evaporate the icy grain mantle. At the 
position of Bl, the outflow shock is the only source that can trigger the ejection of molecules 
from the grain surface. Using the known characteristics of the LI 157 outflow shocks, we 
can deduce which of the two mechanisms is most likely responsible for the formation of the 
detected complex molecules. Multi-transition millimeter observations of NH3, CO and SiO 
indicate that the gas associated with Bl has a temperature between 80 and 300 K and a 
density of about 3 x 10^ cm~^ (Tafalla & Bachiller 1995; Umemoto et al. 1999; Hirano & 
Taniguchi 2001; Nisini et al. 2007). The shocked region associated with Bl extends no more 
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than 20" along the outflow axis (see H2 image in Zhang et al. 2000), and Nisini et al. (2007) 
indicate that shocks in LI 157 have velocities larger than ~ 30 km s~^. We, therefore, 
estimate the time it took the shock to transverse the area associated with Bl to be no more 
than about 1400 yrs. Molecular gas heated by a shock with a velocity of about 30 km s~^ 
takes a few hundred years to cool down below 100 K (Bergin et al. 1998). Hence, the gas in 
the Bl region has been above 100 K for no more than 2000 yrs. 

The Bl region is hot enough for grain mantles to evaporate and it is dense enough for gas- 
phase chemical reactions that result in complex molecules to occur (e.g., Millar et al. 1991). 
However, gas-phase models predict the maximum abundance of complex molecules in the 
gas to occur between a few lO'^ yrs and a few 10^ yrs after the parent molecules are released 
into the gas phase (see Millar et al. 1991) — considerably longer than the timescale estimated 
above — and the maximum abundance of HCOOCH3 predicted by these models is a factor of 
ten less than our observed abundance. Moreover, recent results by Horn et al. (2004) indicate 
that gas-phase production of HCOOCH3 is much less efficient than previously considered by 
gas-phase models. It is highly improbable that most of the observed complex molecules in the 
LI 157 outflow are produced in the gas phase. We therefore conclude that the relatively high 
abundance of complex molecules in the LI 157 outflow is better explained by the formation of 
these species on the grain surface and their subsequent release into the gas phase caused by 
the outflow shock. We note that our results do not discard the possibihty of small differences 
in the formation mechanism of these molecules, and that a small fraction of some of the 
species may exist in the gas phase independent of outflow shocks. For example, HCOOH has 
been observed in a quiescent dark cloud, yet HCOOCH3 has only been observed in active 
regions (Turner et al. 1999; Rcqucna- Torres et al. 2007). The derived HCOOH abundance in 
these studies is about 10^^'^, two orders of magnitude smaller than our abundance estimates 
in LI 157. Hence, only a negligible fraction of the complex species observed in L1157 could 
be due to processes unrelated to the outflow shock. 

Recent surveys of abundances of complex species in hot cores, hot corinos, (the pre- 
sumed low-mass counterparts of hot cores), and molecular clouds in the galactic center (GC) 
region favor the grain surface formation scenario (e.g., Bottinelli et al. 2007; Bisschop et al. 
2007; Requena- Torres et al. 2006; 2007). The abundance ratio of the complex molecules with 
respect to CH3OH (another, more abundant, grain mantle constituent) can be used to inves- 
tigate the origin and evolution of these complex species (e.g., Bottinelli et al. 2007). Using 
the abundance of CH3OH in Bl reported by BP97 (about 10~^), we derive the abundance 
ratios HCOOCH3/CH3OH, HCOOH/CH3OH, C2H5OH/CH3OH, and CH3CN/CH3OH to 
be in the order of 10~^, 10~^, 10^^, and 10~^ respectively. The flrst three are in agreement 
(within an order of magnitude) with the average abundance ratios found in hot cores and 
molecular clouds in the GC region (Requena- Torres et al. 2006; Bisschop et al. 2007), and are 
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also consistent with the upper hmits obtained toward the L1448 outflow by Requena-Torres 
et al. (2007). Taken together these results suggest the dust in hot cores, GC molecular 
clouds, and the LI 157 molecular cloud have similar mantle composition (as also argued by 
Requena-Torres et al. 2006; 2007). In hot corinos, unlike the other sources, the abundance 
of HCOOCH3 and HCOOH with respect to methanol is two orders of magnitude higher that 
in the Bl position of LI 157. Such high relative abundance in hot corinos could be due to 
a longer warm-up phase of the gas surrounding low-mass protostars, when complex species 
are produced relatively more efficiently (Garrod & Herbst 2006). The CH3CN to CH3OH 
ratio is signiflcantly lower by about three orders of magnitude for Bl than for hot cores and 
hot corinos (it has not been measured for GC clouds) . One possible explanation for the low 
CH3CN/CH3OH in Bl could be that processes in the shocked region (but not present in hot 
cores or hot corinos) rapidly destroy CH3CN once it is in the gas phase. It is also possible 
that CH3CN is truly a daughter specie and it shows very low abundance because there has 
not been enough time (< 2 x 10^ yr) for large abundances of this molecule to form in the 
warm gas associated with Bl. With our current data we cannot confidently state which 
scenario is more likely and further observations are needed. 

Similar abundance ratios of a number of molecules in the LI 157 outfiow (and in hot 
cores) compared to the abundance ratio of molecules in the comet Hale-Bopp led Bockelee- 
Morvan et al. (2000) to argue that there is a direct link between cometary and interstellar 
ices. The estimates of HCOOCH3/CH3OH and HCOOH/CH3OH we obtain for L1157 are 
roughly similar to those of Hale-Bopp. Our results are consistent with the conclusions 
reached by Bockelee-Morvan et al. (2000) that molecules in cometary ices could have formed 
by processes very similar to those that produce the chemically rich icy mantles on interstellar 
grains. 

In summary, our results clearly indicate that in molecular clouds with only low-mass 
star formation complex organic molecules can form through grain surface reactions. Outflow 
shocks can heat the surrounding medium and evaporate the icy grain mantles where the 
complex species reside, thereby releasing them into the gas phase and chemically enriching 
the circumstellar environment. Our results show that a protostar's radiation is not the sole 
mechanism that can generate complex molecules near forming stars and that the impact 
of outflows needs to be considered when studying complex species around protostars. If 
no subsequent chemical reactions alter the abundance of complex molecules once they are 
released into the gas phase by the shock, then the abundance estimates from millimeter ob- 
servations can be used to study the mantle composition of the dust in the cloud. Comparing 
the results from our observations of the Bl clump in the LI 157 outflow with those of other 
regions observed by others suggest that the grain mantle composition in the L1157 dark 
cloud is comparable to that of the grains in hot cores and molecular clouds in the Galactic 
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center region. Observations of more complex molecules and estimates of their relative abun- 
dance towards other chemically active outflows will allow us to determine the reliability of 
millimeter observations of shocked molecular gas for estimating the abundance of complex 
molecules in grains and if similar grain mantle compositions are found in different regions of 
low-mass star formation. 
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Table 1. Detected Emission Lines 



Molecule 


Transition 
Line 


Frequency 
[MHz] 


Eu/k 
[K] 


[km s~^] 


rp b 

[mK] 


[km s~^] 


[K km s"^] 


rms 
[mK] 


HCOOCH3-A 


72,5 — 62,4 


90156.48 


19.68 


2.1 


11 


5.0 ±1.5 


0.06 ± 0.02 


2 




80,8 — 7o,7 


90229.63 


20.07 


1.0 


12 


3.7 ± 1.5 


0.05 ±0.02 


3 




83,6 — 73,5 


98611.15 


27.26 


1.0 


15 


4.5 ±0.9 


0.07 ±0.02 


3 




84,5 — 74,4 


98682.60 


31.90 


3.0 


7 


3.3 ±2.5 


0.03 ± 0.03 


2 




202,19 — 102,18 


226718.70 


120.27 


0.8 


26 


2.3 ±0.8 


0.07 ±0.03 


8 


HCOOCH3-E 


72,5 — 62,4 


90145.72 


19.69 


1.0 


12 


5.7 ±1.0 


0.07 ± 0.02 


2 




80,8 — 7o,7 


90227.63 


20.10 


1.0 


12 


3.0 ±0.7 


0.04 ± 0.01 


3 




83,6 — 73,5 


98606.85 


27.28 


1.0 


11 


3.3 ± 1.2 


0.04 ±0.02 


3 




84,5 — 74,4 


98712.06 


31.91 


3.0 


7 


5.5 ±2.3 


0.04 ± 0.02 


2 


CH3CN 


143 - 133 


257482.80 


156.77 


1.2 


26 


4.7 ± 1.2 


0.13 ±0.04 


7 




14i - 13i 


257522.50 


99.89 


1.2 


23 


4.0 ± 1.0 


0.10 ±0.03 


7 




14o - 13o 


257527.40 


92.78 


1.5 


34 


3.6 ±1.0 


0.14 ± 0.04 


6 


HCOOH 


42,2 — 32,1 


90164.63 


23.54 


3.3 


8 


7.3 ±1.8 


0.06 ± 0.02 


2 




62,4 - 52,3 


135737.76 


35.48 


1.4 


17 


5.8 ± 1.4 


0.11 ±0.03 


3 


C2H5OH 


4l,4 - 3o,3 


90117.61 


9.36 


3.3 


10 


6.9 ± 1.0 


0.08 ± 0.02 


2 



^Width of velocity channel. 

^Peak intensity from Gaussian fit to line. 

'^Velocity width and error from Gaussian fit to emission hne. 

'^Errors in the integrated intensity estimate come from propagation of errors, using the error in 
the velocity width and the rms. 
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Table 2. Estimates of Temperature, Column Density and Abundance 



Molecule 


Trot 


N 


X = N/Nh, 




[K] 


[10^^ cm-2] 


[io-«] 


HCOOCHa-A^ 


27 ±4 


15 ±4 


11 ±3 


HCOOCHg-E^ 


18 ±13 


12 ±11 


8±7 


CHsCN^ 


110 ± 50 


0.1 ±0.05 


0.07 ± 0.04 


HCOOH'^ 


10 


8 


5 




25^^ 


10 


7 



^Error estimates come from linear fit to rotational diagram. 

^No error estimate included as only two points were used for linear fit to rotational diagram. 

^No error estimate included as column density was not obtained using the rotational diagram, 
see text. 

'^Trot assumed to be similar to the one derived for HCOOCH3. 
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Fig. 1. — Sample spectrum from our observations of the Bl position in the L1157 outflow. 
The spectrum shown here is centered around 90.17 GHz, spans about 0.15 GHz, has a 
spectral resolution of 0.625 MHz (~ 2 km s~^) and an rms of 2 mK. 
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Fig. 2. — Rotational diagrams of molecules for which we detect two or more transitional 
lines. Errors shown are derived from the errors in the integrated intensity (see Table 1). 



